SPECTRAL THEORY FOR SECOND-ORDER VECTOR EQUATIONS 
ON FINITE TIME- VARYING DOMAINS 

DOUGLAS R. ANDERSON 

Abstract. In this study, we are concerned with spectral problems of second-order vec- 
tor dynamic equations with two-point boundary value conditions and mixed derivatives, 
where the matrix-valued coefficient of the leading term may be singular, and the domain 
is non-uniform but finite. A concept of self-adjointness of the boundary conditions is in- 
troduced. The self-adjointness of the corresponding dynamic operator is discussed on a 
suitable admissible function space, and fundamental spectral results are obtained. The 
dual orthogonality of eigenfunctions is shown in a special case. Extensions to even-order 
Sturm-Liouville dynamic equations, linear Hamiltonian and symplectic nabla systems on 
general time scales are also discussed. 



1. Introduction 

Shi and Chen [27] provide an analysis for second-order vector difference equations of 
the form 

-V (C n Ax n ) + B n x n = \cj n x n , ne[l,N]nZ, N>2, 

where the forward difference operator Ax n := x n +i — x n and the backward difference 
operator Vx n := x n — x n -i are utilized. This is one vector form related to the well-known 
self-adjoint scalar form, see Kelley and Peterson [20, Chapters 6 — 7], 

-A (c n _iAx n _i) + b n x n = Xu n x n . 

Recently, however, there has been a growing interest in time- varying systems such as by 
Davis, Gravagne, Jackson, and Marks |12| or Gravagne, Davis, and DaCunha [14] , that is 
discrete systems with varying step-size between domain points; these systems are sometimes 
also called isolated time scales. An important example in this class of problems would be 
quantum (^-difference) equations; see Simmons [311 Chapter B.5] and Kac and Cheung |_19] 
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for quantum calculus, and Erbe and Hilger |13j and Bohner and Peterson [TO] for more on 
time scales, including isolated time scales. 

Motivated by the appeal of time-varying domains and the ability to simultaneously 
unify and generalize recent and classic results, we introduce here a finite-dimensional anal- 
ysis of second-order vector dynamic equations of the form 

(1.1) - (Px A ) v (t) + Q(t)x(t) = \uj(t)x(t), t £ [a, b]j, b > a(a), 



with boundary conditions 

<"> +s (m^)=°- 

A solution x of (jl.ip is defined on [p(a),a(b)]j. Here T is a (finite) isolated time scale, the 
discrete interval is given by 

[a,b]j := {a, a {a), a 2 (a), ■ ■ ■ ,p(b),b}, b = a N ~ l (a), 

the d X d matrices Pit) for t £ [/o(a), b]i, Q{t) and uj{t) for t € [a, b]j, are all Hermitian with 
P p (a) and P(b) invertible, uj(t) > (positive definite) for t € [a, 6]t, and R and S in (jl.2p 
are 2d x 2d matrices with rank(i?, S) = 2d. The differential operators in (jl.ip are given, 
respectively, by 

xA{ x°(t)-x(t) ^ x v {t) = x{t)-af(t) =xA 

where time-varying step-sizes are given by p a {t) := a(t) — t > 0, n p (t) := t — p{t) > 0, 
with a(t) the immediate domain point to the right of t, p{t) the immediate domain point 
to the left of t, and where compositions are often denoted f(a(t)) = f a (t). Note that the 
coefficient matrix P(t) in (jl.ip may be singular for t € [a,p(6)]ir. Moreover, the second- 
order dynamic operator given in (jl.ip may not be formally self-adjoint in the space of vector 
series of the form 

{^»)}n=-i, v-Ha) := p(a), a°(a):=a, a n {a) := a{o n ~\a)), 



subject to the boundary conditions in (jl.2p unless those conditions themselves are self- 
adjoint in some sense to be defined. On the other hand, even if the boundary conditions 
(jl.2p are self-adjoint, the second-order dynamic operator may still not be self-adjoint. As 
in the difference equations case [27], self-adjointness depends on whether the vectors at the 
boundaries, namely x p {a) and x a {b), can be solved from (jl.2p in terms of x{a) and x(b). 
Indeed if this is the case, we will call (jl.2p proper self-adjoint boundary conditions. If 
(|1.2p happen to be improper, we can find a smaller space where the second-order dynamic 
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operator is self-adjoint. Once this is sorted out, we will give the appropriate analysis of the 
eigenvalue problem for (11. lj) and (11. 2p . 

The particular appeal of (jl.ip is that it is still a discrete problem, but with non-constant 
step-size between domain points. As 

lim x A (t) = lim x v (t) = x'(t), ' = — , 
/v^o /v->o at 

for differentiable functions x, these dynamic results serve as an alternative discrete analog to 
the differential equations case. There are many papers on the spectral theory for difference 
equations, including recent papers such as Ji and Yang [T71 [18], Shi [25], [26], Shi and Chen 
[27] [28], Shi and Wu [29], Shi and Yan [30], Sun and Shi [32], Sun, Shi, and Chen [33], Wang 
and Shi [Mj, and Bohner, Dosly, and Kratz [SJ, but none on dynamic vector equations. Thus 
this work will continue and extend the discussion of self-adjoint equations on time scales 
found in [¥[ [5j \TT \ \15 \ I23|. As in the uniformly-discrete case [27], we use mixed derivatives, 
as —V is a natural adjoint of A. To emphasis the dimensional analysis we focus on isolated 
time scales, but a similar treatment can be made on general time scales. 

The analysis of (jl.ip and its solutions will unfold as follows, largely motivated by Shi 
and Chen [27] . In Section 2 we will give a definition of the self-adjointness of the boundary 
conditions in (|1.2|) . Section 3 contains the introduction of a suitable admissible function 
space in which the corresponding dynamic operator is self-adjoint, followed in Section 4 by 
some fundamental spectral results. In the event of proper self-adjoint boundary conditions, 
the dual orthogonality of eigenfunctions will be given. In Section 5 we discuss the possibility 
of extending these results to even-order Sturm-Liouville equations, linear Hamiltonian and 
symplectic nabla systems on time scales. 

2. Self-Adjoint Boundary Conditions 
In this section we are concerned with the second-order vector dynamic operator I given 

by 

(2.1) £x(t) = uo-\t) [-(Px A f(t) + Q(t)x(t)] , t € [a, %, 

where we will write x E S% if {x(a n (a))}^ = _ 1 satisfies the boundary conditions given in 
(fL2D . If we set 

e[p(a),a(b)} := {x = {x(a»)}£U : x(t) £ C d , t G [p(a), a(b)])r} , 
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then dim £[p(a), a (b)] = d(N + 2). Define the (weighted) inner product to be given by 
(2.2) 



\X,v) = ^2 y*( t ) UJ {t)x{t)pp{t) for x,y G £[p(a),a(b)], 

te[a,b] T 



where u(t) is as in (jl.ip . p p (t) = t — p is the left-graininess at t, and y*(t) denotes the 
complex conjugate transpose of y(t). Using standard notation, x _L y will mean (x,y) = 0. 

Theorem 2.1 (Lagrange Identity). Assume the dxd matrices P(t) fort € [p(a),b]j, Q(t) 
and u>(t) for t £ [a, b]r, are all Hermitian with P p {a) and P(b) invertible, and co(t) > for 
t £ [a, b]j. Then for x, y € £[p(a), cr(b)] we have that 

(2.3) (£x,y) - (x,£y) = [(Py A )*x - y*Px A ] (b) - [(Py A )*x - y*Px A ] P (a), 



where the inner product is given in (12. 2p . 

Proof. Let x,y € £[p(a),a(b)]. Using summation by parts [U Theorem 2.7] we have 
(£x,y) - (x,£y) = ^ y*(t)u;(t)£x(t)p p (t) - ^ (£y)* (t)u(t)x(t)p p (t) 

t£[a,b] T te[a,fc] T 

= - E y*(t){(Px A f-Qx)(t)p p (t) 

te[a,b]t 

+ £ ((Py A ) v -Qy)*(tM% p (t) 

te[a,% 

= - £ y*(t)(Pz A ) v (t) AV (t) + ^ (Py A ) v *(t)x(t)/x p (t) 

te[a,&] T tG[a,fc] T 

= - [y*Px A ] (t) b + J2 {Py A T P {t)x V (t)p p (t) 
PW t 5~M 

+ (Py A f*(t)x(t)p p (t) 

te[a,b]t 



[y*Px A ] (t) + [(Py A yx] (t) 
PW 



p(a) 



[{Py A Tx-y*Px A ] (t) 



b 

P (a) 



and the result follows. 



□ 



The above theorem, Theorem 12.14 is a ls° valid on general time scales; see Anderson 
and Buchholz [UJ Theorem 2.13]. 



Definition 2.2. The boundary conditions given in (jl.2p are called self-adjoint iff 

b 



(2.4) 



[(Py A )*x - y*Px A ] (t) 



p(a) 
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for x, y G 8%. 

The next theorem then follows easily from the definition and Theorem 12,11 

Theorem 2.3. Assume the d x d matrices Pit) for t G [p{a),b\j, Q(t) and uj(t) for t G 
[a, b]j, are all Hermitian with P p (a) and P(b) invertible, cu(t) > for t G [a, b]j, and R and 
S in (|1.2|l are 2d x 2d matrices with rank(i?, S) = 2d. If the boundary conditions in (jl.2p 
are self-adjoint, then 

(£x,y) = (x,£y) 

for £[p(a),a(b)] with x, y G 8%. 

Lemma 2.4. Assume that R and S are 2d x 2d matrices with rank(i£, S) = 2d. Then 

(2.5) SR* = RS* 

if and only if the boundary conditions (jl.2p are self-adjoint. 

Proof. We follow Kratz [221 Proposition 2.1.1]. Let x G £[p(a),a(b)] with x G M. By the 
definition of x satisfies the boundary conditions (|1.2p . Let /C = (^) be any 4d x 2d 
matrix with Im/C = ker(l?,C), so that SB + i?C = 0, rank/C = 2d, and so that x G 8$ 

if and only if C^(b)j = -^ r /' (^(6)^(6^) = ^ ^ or some *7 e ^ 2d - This results in (|2.4p 
holding for all x, y G 8? if and only if 77* [£>*C — C*B]r]2 = for all 771, 772 G C 2d , that is to 
say B*C = C*B. 

First, assume that RS* = SR*. Since rank(i?, S) = 2d, we see that Im ( ) = 
ker(i2, S). Consequently, the matrix K, given above can be taken to be ( J* R * ) ; by the Her- 
mitian assumption on RS* we see that B*C is Hermitian, and thus the boundary conditions 
(|1.2p are self-adjoint by Definition 12.21 The rest of the proof is identical to |22} Proposition 
2.1.1] and is omitted. □ 

Lemma 2.5. Assume the dxd matrices P(t) fort G [p(a),b]j, Q(t) andu{t) fort G [a, b]i, 
are all Hermitian with P p (a) and P(b) invertible, u)(t) > for t G [a, b]j, and R and S in 
(jl.2p are 2d x 2d matrices with rank(i?, S) = 2d. If the boundary conditions in (jl.2p are 
self-adjoint, then x G 8$ if and only if there exists a unique vector rj G C 2d such that 

Proof. See Kratz [22\ Proposition 2.1.2] and Shi and Chen |27l Lemma 2.2]. □ 
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3. Admissible Function Space and Self-Adjointness of the Dynamic Operator 



Throughout this section we assume the d x d matrices P(t) for t G \p(o), b]i, Qif) an d 
ui(t) for t G [a, b]j, are all Hermitian with P p {a) and P(b) invertible, u(t) > for t G [a, b]j, 
the matrices R and S in (jl.2p are 2(i x 2d matrices with rank(i?, S) = 2d, and the boundary 
conditions in (jl.2p are self-adjoint. We begin the construction of an admissible function 
space to establish the self-adjointness of a dynamic operator related to (jl.ip . following the 
development in the uniformly discrete case presented in [27J. 

Let R = (Ri, R2) and S = (Si, S2), where Rk and Sk are 2d x d matrices, k = 1,2. By 
Lemma 12.51 x G M if and only if there exists a unique vector i] G C 2d such that 

x p {a) = Stv, x(a) = {St + /x p (a)P"- 1 (a)^)r ? , 
(3.1) x(b) = -S* 2V , x°(b) = (-S* 2 + ^(b)p-\b)R^) V , 

where P p ~ l = Rewrite the boundary conditions (jl.2p as 

where the coefficient matrix V is given by 

r = (rxPO- 1 ^) + -4t5i, -^tS 2 ) . 

V /z p (a) /i CT (6) / 

Set 

r = rank f fl x + -J—S^^a), ^77^2) G {0, 1, • • • , 2d}, 

where this value r will play a key role in the development to follow. If the matrix V 
is invertible, in other words if r = 2d, then the boundary conditions f)1.2|) are proper, 
since (|3.2p is solvable for x p (a) and x a (b), respectively. We then see that, as x p (a) and 
x a (b) are not weighted with respect to the weight function 00, the diV-dimensional space 
L[p(a),a(b)] := {x G i[p(a),a(b)\ : x G might serve as a suitable admissible space. If 
r < 2d, however, then from (|3,2p we have that the 2d components of x(a) and x(b) are 
themselves linked by 2d — r relations. From (jl.ip we see that x(a) and x(b) are weighted 
by uj(a) and u:(b) in the first and last vector equations. In the two vector equations, via 
some transformation, only r scalar equations are really weighted, while the remaining 2d — r 
ones involve x p (a), x(a), x a (a), x p (b), x(b), and x a (b) but not A; consequently they can be 
viewed as extra conditions for the admissible functions. We will now show this as follows. 
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Using standard matrix theory [16] , there exist 2d x 2d-unitary matrices U and V such 

that 



(3.3) 



U*FV = diag{0,M}, 



where M is an r x r matrix with rankM = r. Let U = (U\, Uq) and V = (V\, V 2 ), where 
U\ and V\ are 2d x {2d — r) matrices, U 2 and V 2 are 2d x r matrices. Since V is unitary, 
from (13.31) we have that 



(3.4) 



(3.5) 



V{Vi=I 2 d-r, V{V2 = 0(2d-r)xr> ^2*^2 = Ir, 



U-lF — 0(2d-r)x2d 



, V 2 = F*U 2 M 



*-i 



Using (|3.ip and (|3.3p . we see that 

(3-6) = Udiag{0,Ar}t/*7 ? = (0 2dx(2(i _ r) ,U 2 M* )C/*r ?J 

so that by (|3.4|) we have 



(3.7) 



Vj* (x T (a), — x T (b)) = 0, x T = transpose of 



x. 



The first and last relations in (jl.ip can be written as 
( \ ( \ ( 



A 



x(a) 



-x(b) 
V 7 



oj- l (a] 



V 



/ 



u- l (b) 
where we have taken 

P(t) = Q(t) + 



^PP(a)xP(a) + P(a)x(a) - 1 ^^P(a)x°(a) 



From (|3,7p we obtain the 2d — r scalar equations 
(3.8) 



( 



v* 



\ 



uj- l (a) 



V 



u- l {b) 



( 



(a 



^ PP (a)xP(a) + P{a)x(a) - ^P(«K(a) 



V 



^P^(b) - P(b)x(b) + 1 ^^ ) P{b)x°{b) 



If x satisfies (|3.8[) we will write i£ / Noting the absence of A from (|3.8p . we can think of 
(|3,8p as additional conditions together with the boundary conditions (jl.2p . Thus, we define 
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the admissible function space 

L l[p(a), <r(b)] := {x G £[p(a),a(b)] :i£fn4 



From (|3,4p we see that x G stf if and only if x satisfies 



(3.9) 



\ 



-i 



Op(a)) 



^(a)xP(a) + P(a)x(a) - -^-^ P ( a )^(a 



^POib^ib) - P(b)x(b) + j^^P^Cft) 



\ 



w(o) 
u(b) 



y 2l 



for some 7 G C r . As P p {a) and P(b) are invertible, (j3.9[) implies that 



2^ o 



V 7 



\ 



PP-\a) U» p (a)y P(a)x(a) - gjgP(a)s»(a) 



(-gjgP^6)x"(6) +/i p (6K(6)P(6)x(6) 
/ 



(3.10) 



+ 



( Mp (a)) 2 P"- 1 (a)a;(a) 



V 







V 2 j. 



Let x G L 2 [p(a), cr(b)]. Inserting (|3.10p into boundary conditions (jl.2p or (|3.2p yields 

/ ^ 



r diag{(/i p (a)) 2 w(a), // p (6)/x CT (6)a;(6)}V r 27 

(3.11) 
(3.12) 
From (13.31) we have 



W«)) 2 P(«)x(fl)-jgp(«)/(« 



^^(6)^(6) - n p {b)n a (b)P{b)x(b) 



■ 1 S 1 P p (a),R 2 -^-S 2 P(b)) (» 



Mp( 

=:/(x(a),x CT (a),x"(6),x(6)). 



x(a) 
(6) 



(3.13) 



rdiag{(^ p (a))" ! w(a),/i p (6)/i (T (6)a;(6)}V2 = U (0 rx(2d _ r ), J T ) T , 



where J = MV 2 * diag{(^ p (a)) a; (a), /x p (6)//o-(6)w(6)}V2 is an r x r invertible matrix, as 
rankM = rank V 2 = r , and a; (a) and oj(b) are positive definite by assumption. Multiplying 
(|3.12p on the left by U* , we have from (|3.5p (in particular f7*fi»2 = 0) and (|3. 13[) that 



(3.14) 



H p {a 



■SiP p {a)x{a) + R 2 x(b) = 



(3.15) 
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7 = J~ X Vtf (x(a), x», x»(b), x(b)) . 



We are now in a position to get a useful characterization for a function x to be a member of 
L^[p(a),a(b)]. As we have the necessary structure in place and in parallel with the discrete 
case presented by |27j . we omit the proofs of the following key results; for more details, 
please see [27] , 

Lemma 3.1. We have x € L^[p(a),a(b)] if and only if x satisfies (|3.1Up and (|3.14p in 
which 7 is determined by (|3.15|) , 

Lemma 3.2. We have m := dim L^[p(a) , cr(b)} = d(N — 2) + r. 

Theorem 3.3. Assume the d x d matrices Pit) for t € [p(a),b]x, Q(t) and uj(t) for t € 
[a, 6]tj o,re all Hermitian with P p (a) and P{b) invertible, oo{t) > for t € [a, 6]t, and R 
and S in (jl.2p are 2d x 2d matrices with rank(i?, S) = 2d. If the boundary conditions in 
(jl.2p are self-adjoint, then L^[p{a),a(b)] is an m-dimensional Hilbert space with the inner 
product defined in (|2,2p . 



Let 



L := 



\Ll[p(a),a{b)} 



, where i is given in (|2.ip . 



We will see the self-adjointness of this dynamic operator. 

Lemma 3.4. The dynamic operator L maps L^[p(a),a(b)] into itself. 

Proof. See [271 Proposition 3.3], which uses (|3.10p with x replaced by Lx as in 



/ \ 

LxP(a 



Lx a (b) 
V 7 



(3.16) 



P^(a) Up p (a)f P(a)Lx(a) - %$P(a)Lx° (a) 



{~lfM PP{b)LxP{h) + ^)^(b)P{b)Lx(b) 



( 



+ 



\ 



ipi p {a)) 2 P^iaMa) 



\ 







/ u p (6)^ (T (6)P~ 1 (^(6) 



V27, 



where 7 is determined by (|3.15p with x(a), x a (a), x p (b), and x(b) replaced by Lx(a), Lx a {a), 
LxP(b), and Lx(b). □ 
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Theorem 3.5. Assume the d x d matrices Pit) for t E [p{a),b\j, Q(t) and uj(t) for t € 
[a, b]j, are all Hermitian with P p (a) and P(b) invertible, u)(t) > for t £ [a, b]j, and R and 
S in (jl.2p are 2d x 2d matrices with rank(i?, S) = 2d. If the boundary conditions in f)1.2|) 
are self-adjoint, then the dynamic operator L is self-adjoint. 

4. Dynamic Spectral Theory 

As in the previous section, the machinery for a dynamic spectral theory is in place. All 
proofs of the following results are modifications of those found in the discrete case in |27| . 
and thus are omitted here. For completeness, however, we list the main ideas and theorems 
that hold in this setting as well. 

Definition 4.1. A complex number A E C is an eigenvalue of (|l.ip . f)1.2|) if and only if 
there exists x E £[p(a),a(b)] that is nonzero and solves the boundary value problem (|l.ip . 
(|l,2p . The nonzero solution x is a corresponding eigenfunction for X, denoted x(-,X). 

Following the discussion in the previous section, any eigenfunction of (jl.ip . (jl.2p is 
in the space L"^[p(a), o~(b)]. Therefore, the following fundamental spectral results for (jl.ll) . 
(jl.2p can be directly concluded by Theorem 13.51 and by employing the spectral theory of 
self-adjoint linear operators in Hilbert spaces, spaces. 

Theorem 4.2. Assume the d x d matrices Pit) for t G [p(a),b]j, Q(t) and u(t) for t G 
[a, b]x, are all Hermitian with P p (a) and P(b) invertible, u(t) > for t € [a, b]x, R and S 
in (|1.2p are 2d x 2d matrices with rank(ii, S) = 2d, and the boundary conditions in (|1.2p 
are self-adjoint. Let r = rank ^Ri + - ^ Si P p (a), j^^S^j and m = d(N — 2) + r. 

(i) The eigenvalue problem (|l.ip . (II. 2p has m real eigenvalues 

(4.1) Ai, A2, • • • , A m , 

including multiplicities, and m linearly independent eigenf unctions 

(4.2) x(-,Ai),x(-,A 2 ),--- ,x(;Xm), 
that form an orthonormal set, to wit 

(4.3) (x(; Xj),x(; Ajfc)) = ^2 ^{t^k^i^xit^X^pvit) = 5 jk , j,k G {1,2, ••• ,m}. 

te[o,6] T 
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(ii) The eigenfunction basis for (jl.ip . (jl.2p consists of m linearly independent eigen- 
functions (14, 2p and is complete for admissible function space L^[p(a),a(b)], to wit 
for each x G L^[p(a), o~(b)], there exists a unique set of scalars {cj}JL 1 C C such 
that, for t G [p(a),a(b)]j, we have 



(4.4) x(t) = ^2cjx(t,Xj), 

(4.5) Cj = (x,x(-, A,-) = ^ x*(t,\j)u}(t)x(t)n a (t), j G {1,2, • • • ,m}, 

te[o,6] T 

and i/ie following Parseval equality holds, namely 

m 

(4.6) <x,x) = J>/. 

i=i 

(iii) T/ie dynamic operator L has the spectral resolution 

771 roc 

(4.7) Lx(t) = J2 X W x ( t )= XdE x x(t), xeLl[p{a),a(b)], t G [p{a), <r(b)] T , 

3=1 J -°° 



(4.8) 



with Lx p (a) and Lx a (b) given as in (|3.16p . the projective operators are given by 

Kjx{t, Xj) = 5 jk x(t, A fe ), j, fc G {1, 2, • • • , m}, t e [p(a),a(b)]j, 
and the projective operator-valued function is given by 



(4.9) 



Ex 



zJa<A,<i> 



: A > 0, 
7T,-, : A < 0. 



Theorem 4.3 (Dual Orthogonality). Assume the d x d matrices P(t) for t G [p(a),b]j, 
Q(t) and u)(t) for t G [a, 6]t> «re all Hermitian with P p (a) and P{b) invertible, ui(t) > 
for t G [a, b]x, R and S in (jl.2p are 2d x 2d matrices with rank(i?, S) = 2d, and the 
boundary conditions in (jl.2p are self-adjoint. Let r = rank (^R\ + - j^ S\P p (a), - 
and m = d(N — 2) + r. 7/ i/te boundary conditions (|1.2p are proper, that is r = 2d, then 
m = dN and the dN eigenf unctions f|4.2[) in Theorem \4-%\ satisfy the dual orthogonality 
condition 



(IN 



(4.10) J2x(t,X j )x*{r,X j ) = 6{t,T)uj- 1 (T), t,rG[a,6] T , S(t,r) = < 

3=1 



:t^r, 

1 : t = t. 
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Remark 4.4. There are additional variational properties of eigenvalues and comparison 
results, including Rayleigh's principle and a minimax theorem, that carry over verbatim 
from the discrete case. We refer the interested reader to \17\ Section 5]. One could also use 
the recent results by Shi and Wu [29J in the discrete case to improve our results above. 

5. Further Extensions to Sturmian Theory 

For the discussion that follows, assume the time scale T is general, that is to say not 
necessarily isolated. In the classic continuous case (T = R), scalar, vector, and matrix 
equations of the form 

-(px')'(t) + q(t)x(t) = 

are embedded in a robust theory that extends to even-order Sturm-Liouville equations, 
linear Hamiltonian and symplectic systems; see, for example, Reid |24j and Kratz |22j . 
It is common on general time scales to consider these same connections using the basic 
second-order scalar equation 

(5.1) -(px A ) A (t)+q(t)x°(t)=0. 

Erbe and Hilger |13j began a Sturmian study of this equation, followed by Agarwal, Bohner, 
and Wong [lj. Recently Kong [21] extended some of these results. To the best of our 
knowledge, however, it is an open problem on how to extend (|5.ip to even-order self-adjoint 
Sturm-Liouville problems using only delta derivatives. 

As seen above there is a natural alternative theory that builds on the scalar equation 

-(px A f(t)+q(t)x(t) =0. 

After Atici and Guseinov introduced the nabla derivative [8] as a left-hand counterpart on 
time scales to the delta derivative, there has been an attempt to formulate a Sturmian 
theory using both delta and nabla derivatives. Messer [23] continued the work of [8] by 
studying second-order scalar self-adjoint equations; this work was extended to higher order 
equations, first by Guseinov [15], and then by Anderson, Guseinov, and Hoffacker [5]; see 
also Davidson and Rynne [11]. A Reid roundabout theorem and an analysis of dominant 
and recessive solutions is given in [61 [7] for the related matrix equation (|5,2p . given below. 

A next step might be to extend the results of this paper to even-order Sturm-Liouville 
equations and linear Hamiltonian systems on arbitrary time scales, or more generally to 
symplectic systems on time scales. We illustrate this in one possible treatment not found 
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in the literature. For example, we reconsider (|2.ip rewritten in the related homogeneous 
second-order self-adjoint matrix dynamic equation in the form 



(5.2) 



PX A f (t) + Q(t)X(t) = 0, te[a,b] T . 



Let P be invertible Hermitian, Q ld-continuous Hermitian, and let B denote the set of all 
n x n matrix- valued functions X defined on T such that X A is continuous on [p(a), b]j and 
(PX A ) V is left-dense continuous on [a, 6]t- Then X is a solution of (|5.2p on [p(a),a(b)]j 
provided IgD and X{t) satisfies the equation for all t € [a, b]j. 



o 



Lemma 5.1. Let X € B. If we set 

(5.3) Z:=( p ^ A ) and S :— , () 
then X solves (j5.2[) on [a, b]j if and only if Z solves 

(5.4) Z^ = S(t)Z, te[a,a(b)] T . 

Proof. We give an idea of the proof. Let X € B and Z, S be given as in (|5.3p ; multiplying 
out, we get 



SZ 



-tipPP^Q pp- 1 \ ( x 



o 



)(pW = ( 



~ f j, p PP- 1 QX+PP- 1 PX l 
QX 



-fj, p PP- 1 (PX A ) v +PP- 1 PX A 
QX 

-PP- 1 ((PX A )-(PX A )P)+PP- 1 PX A 
QX 



PP- 1 (PX A )P 
QX 



\QX ) 



A v 

AW 



{PX A ) 



where (PX^) P = P P X^ since X € B implies X A is continuous. 



□ 



Next we introduce a new type of Hamiltonian for the nabla differential operator, not 
previously discussed in the literature, although clearly in parallel with that available for 
the delta differential operator. In the discrete case see Ahlbrandt and Peterson [3j, and Shi 
, 126] . and see Bohner and Peterson [1CH Chapter 7] for the delta operator on time scales. 



Definition 5.2. A 2n x 2n-matrix-valued function H is Hamiltonian with respect to T 
provided 



(5.5) is Hamiltonian and I + n p (t)rl(t)A4*A4 is invertible for all t 6 T, 
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where we have used the 2n x 2n-matrix M = ((jo)- ^or % satisfying (|5.5j) . i/ie system 

(5.6) z v = ^(t)[X*Mz' , + 7W>I*z] 
is i/ien called a linear Hamiltonian dynamic nabla system. 

Remark 5.3. Take J to be the 2n x 2n-matrix J = q) anc ^ ^ = ( C ) ^ or some 
n x n-matrix-valued functions A, B, C such that B and C are Hermitian. Then Ti* J '+ JH = 
0, so that "H is Hamiltonian, and 

I + li P UM*M={ I QI ^). 

Therefore, condition (|5,5p may be rewritten as % = ( ^ _® » ) , where I — fi p A* is invertible 
and B,C are Hermitian. If we let z = we can rewrite the system (j5.6j) as 

(5.7) x v = A(t)x + B(t)u p , u 7 = C(t)x - A*(t)u p . 

Also note that equation ()5.2[) is equivalent to a Hamiltonian system by taking A = 0, 
B = pp^ 1 , C = Q , x = X , and u = PX^. Finally, note that (|5.7p is the form of the 
discrete Hamiltonian introduced by Ahlbrandt [2]. 

Theorem 5.4. // we set 

(5.8) S = {I + n p HM*Mr 1 H, 

then every Hamiltonian system (|5.6p is also a symplectic nabla system z v = S(t)z. 

Proof We assume first that z solves (|5.6p . Using the simple formula z p = z — fipZ^ 7 and 
observing that M*M + MM* = I, we obtain (/ + ^ p HM*M)z v = Hz, so that §5M 
follows from (|5.5|) . To verify that S defined by (|5.8p is symplectic with respect to T, 
videlicet satisfies condition S*(t)J + JS(t) = n p (t)S*(t)JS(t) for all t G T, note that 

<? = ( ■A-n P B(l-fL P A*)- 1 C Bil-fipA*)' 1 \ 
It is straightforward to check that S*J + JS = [i p S*JS holds on T. □ 

As a final important illustration, consider even-order Sturm-Liouville dynamic equa- 
tions; see a related discussion by Shi and Chen [28] on higher-order discrete Sturm-Liouville 
problems. This is a new formulation not previously presented, but clearly in the spirit of 
the classical Sturm-Liouville analysis in the continuous case; see Kratz )22[ (6.2.1)] and Shi 
and Chen [281 P- 8]. 
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(5.9) 



k=0 



Pn-kU 



(t) 



-If M 



A n_1 V 



(*) + 



P3V 



A 2 V 



V 2 A 



+ (P2y VA )^(t)-(pw A f(t)+ P0 (t)y(t), 



which is formally self-adjoint [5] 7 p n ^ 0. We will show (I5.9P can 6e written in the form of 
(|5.7p . where 



1 : i/j = i + 1, 1 < z < n — 1, 
: otherwise, 



A = (a>ij)i<i,j<n with a,ij = < 

£ = diag jo, ...,0, -p 1 , C = diag{p ,Pi,P2) • • • >p£-i} • 

I Pn ) 

To do this, we introduce the pseudo-derivatives of the function y given by 

.[*] _ „.V fc n ^ l„ ^ „ i ,.[01 _ „.V° 



(5.10) 



If 
y ln] 

.[n+k] 



V n_1 A 

PnV 

Pn- k y vn - k - lA - (y [n+k - 1] ) A 

k 

E(-» 



i=0 



^ ( M V "~ ! ~ 1A 



1 < fc < n - 1, 



(5.11) 



,[2n] 



poy- [y 



[2n-l] 



My. 



From (|5.1U|) and (|5.1ip we have 



(yM) V = y[ fc+1 ], 0<A;<n-2, 



.[n-l] 



1 < A; < n- 1, 



[2 "~ 1]N ) = P0 y-My. 



For more details on the establishment of related (but different) types of formulas, please 
see [H Section 3] . Then using the substitution 



x 





, u = 


/ y p™-l]\ 

ypn-2] 1 






\ yM / 
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and the matrices A, B, and C above, we have that 

x v = A(t)x + B(t)u p , u 7 = C(t)x - A*(t)u p , 

and the example is complete. 

Consider extending the discrete linear Hamiltonian analysis given by Shi [25] to an 
eigenvalue problem for the Hamiltonian nabla system (|5.7p on general time scales (including 
non- isolated time scales), namely 

x w = A{t)x + B(t)u p , u 7 = [C{t)-Xoj(t)]x-A*{t)u p , te[a,b] r 

with boundary conditions 

(5-12) R(S ) )+S(^)=0. 



x(b) J \ u(b) 

Here A, B, C and u are d x d matrices, B and C are Hermitian, ui > is positive definite, 
A* denotes the complex conjugate transpose of A, and R and S are as in the first part of 
the paper, to wit 2d x 2d matrices with rank(i?, S) = 2d. Write 



Ait) B(t) 7' " w - \ o o ; ' " - \i o 

Then H and W are 2d x 2d Hermitian matrices with W(t) > for t £ [a, b]i, and the system 
(|5.7p . (|5.12p can be recast in the form 



(5-13) ^(5g)=W) + AW(t)](^| 

Define the natural dynamic nabla differential operator for (|5,7p . (|5.12p via 
(5-14) e(x,u)(t) = j(^)-H(t)(: p %). 

Clearly this includes (jl.ip . ()1.2p as a special case, if we replace b by a(b) as the right-hand 
endpoint to accommodate u = Px^. As before we introduce the linear space 

£[p(a),b] := {{(x(t),u(t))} t€[p{a)Mj : x(t),u(t) G C d , t € [p(o),6]) T } , 

and the weighted inner product 

(x,y)= f y*{t)u(t)x{t)Vt, (x,u), (y,v) € £[p(a),b). 

Jp(a) 

Then we have the following key result. 

Theorem 5.6 (Lagrange Identity). For all (x,u), (y,v) G £[p(a),b], we have 

c(t) 



b {(y*,vne( x , u ) - e( y ,v)* c P )} (t)vt = (y*(t),v*(t)) j (f )N 

p(a) 



p(a) 
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Proof. Suppressing the variable t, we have 



{ y *y*)t(?,u) 



y*u v + y*Cx - y*A*u p + v p *x v - v p *Ax - v p *Bu p 



v v *x + y*Cx - v p *Ax + y v *u p - y*A*u p - v p *Bu p 



so that when we subtract the second from the first, we obtain 



(y*,v p *)l(x,u)-£(y,vy(* P ) 



y u +v H x • +v x — y u f 



(v*xf - (y*uf. 



The result follows from the fundamental theorem of calculus. 



□ 



A spectral theory for dynamic linear Hamiltonian nabla systems should now be possible; 
the interested reader is encouraged to see |25j for the details in the discrete case. 
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